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Abstract 

The main purpose of this article is to provide a guide to theorems on global properties of 
solutions to the Einstein- Vlasov system. This system couples Einstein's equations to a kinetic 
matter model. Kinetic theory has been an important field of research during several decades in 
which the main focus has been on non-relativistic and special relativistic physics, i.e., to model 
the dynamics of neutral gases, plasmas, and Newtonian self-gravitating systems. In 1990, 
Rendall and Rein initiated a mathematical study of the Einstein- Vlasov system. Since then 
many theorems on global properties of solutions to this system have been established. This 
paper gives introductions to kinetic theory in non-curved spacetimes and then the Einstein- 
Vlasov system is introduced. We believe that a good understanding of kinetic theory in 
non-curved spacetimes is fundamental to a good comprehension of kinetic theory in general 
relativity. 




Update (27 May 2011) 

This is a revised and updated version of the article from 2005. A number of new sections on the 

Einstein- Vlasov system have been added, e.g., on the formation of black holes and trapped surfaces, 
on self-similar solutions, on the structure of static solutions, on Buchdahl type inequalities, on the 
stability of cosmological solutions, and on axisymmetric solutions. Some of the previous sections 
have been significantly extended. The number of references has increased from 121 to 197. 
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1 Introduction to Kinetic Theory 



In general relativity, kinetic theory has been used relatively sparsely to model phenomenological 
matter in comparison to fluid models, although interest has increased in recent years. From a 
mathematical point of view there are fundamental advantages to using a kinetic description. In 
non-curved spacetimes kinetic theory has been studied intensively as a mathematical subject during 
several decades, and it has also played an important role from an engineering point of view. 

The main purpose of this review paper is to discuss mathematical results for the Einstein- Vlasov 
system. However, in the first part of this introduction, we review kinetic theory in non-curved 
spacetimes and focus on the special-relativistic case, although some results in the non-relativistic 
case will also be mentioned. The reason that we focus on the relativistic case is not only that it is 
more closely related to the main theme in this review, but also that the literature on relativistic 
kinetic theory is very sparse in comparison to the non-relativistic case, in particular concerning 
the relativistic and non-relativistic Boltzmann equation. We believe that a good understanding of 
kinetic theory in non-curved spacetimes is fundamental to good comprehension of kinetic theory 
in general relativity. Moreover, it is often the case that mathematical methods used to treat 
the Einstein- Vlasov system are carried over from methods developed in the special relativistic or 
non-relativistic case. 

The purpose of kinetic theory is to model the time evolution of a collection of particles. The 
particles may be entirely different objects depending on the physical situation. For instance, 
the particles are atoms and molecules in a neutral gas or electrons and ions in a plasma. In 
astrophysics the particles are stars, galaxies or even clusters of galaxies. Mathematical models of 
particle systems are most frequently described by kinetic or fluid equations. A characteristic feature 
of kinetic theory is that its models are statistical and the particle systems are described by density 
functions / — f{t,x,p), which represent the density of particles with given spacetime position 
{t,x) e M X E'^ and momentum p € K'^. A density function contains a wealth of information, and 
macroscopic quantities are easily calculated from this function. In a fluid model the quantities that 
describe the system do not depend on the momentum p but only on the spacetime point {t,x). A 
choice of model is usually made with regard to the physical properties of interest for the system 
or with regard to numerical considerations. It should be mentioned that a too naive fluid model 
may give rise to shell-crossing singularities, which are unphysical. In a kinetic description such 
phenomena are ruled out. 

The time evolution of the system is determined by the interactions between the particles, which 
depend on the physical situation. For instance, the driving mechanism for the time evolution 
of a neutral gas is the collision between particles (the Boltzmann equation). For a plasma the 
interaction is through the electromagnetic field produced by the charges (the Vlasov-Maxwell 
system), and in astrophysics the interaction is gravitational (the Vlasov-Poisson system and the 
Einstein-Vlasov system). Of course, combinations of interaction processes are also considered but 
in many situations one of them is strongly dominating and the weaker processes are neglected. 

1.1 The relativistic Boltzmann equation 

Consider a collection of neutral particles in Minkowski spacetime. Let the signature of the metric 
be (— , -t-, -f ). In this section we assume that all the particles have rest mass m = 1, and we 
normalize the speed of light c to one. We point out that in Section [2] on the Einstein-Vlasov 
system, the dependence on the rest mass and the speed of light will be included in the formulation 
of the system. The four-momentum of a particle is denoted by p", a = 0, 1, 2, 3. Since all particles 
have equal rest mass, the four-momentum for each particle is restricted to the mass shell, p°'Pa = 
—rn^ = — 1. Thus, by denoting the three- momentum hy p E R^, p'^ may be written p"^ — {p'^,p), 
where p^ = -^/l + |pp is the energy of a particle with three-momentum p, and |p| is the usual 
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Euclidean length of p. The relativistic velocity of a particle with momentum p is denoted by p and 
is given by 

P= J-^ - (1) 

Note that |p| < 1 = c. The relativistic Boltzmann equation models the spacetime behavior of the 
one-particle distribution function / = f{t,x,p), and it has the form 

dt + ^-^?jf^Q{fJ), (2) 

where the relativistic collision operator Q{f,g) is defined by 

Qif,9)^ / k{p,q,uj)[f{p + a{p,q,uj)uj)g{q-a{p,q,uj)uj)- f{p)g{q)]dujdp. (3) 

Note that g = f in Equation Here duj is the element of surface area on §^ and k{p, q, uj) is the 
scattering kernel, which depends on the differential cross-section in the interaction process. We 
refer to |178] . [53] and [55] for examples of differential cross-sections in the relativistic case. The 
function a{p,q,u}) results from the collision mechanics. If two particles, with momentum p and q 
respectively, collide elastically with scattering angle w G S^, their momenta will change, i.e., p ^ p' 
and q ^ q' . The relation between p, q and p', q' is given by 

p' = p + a{p,q,uj)uj, q' = q- a{p,q,uj)uj, (4) 

where 

2(p" + g>V(c..(g-p)) 
"(^''^'"^^ ipO+qO^-i.-ip + qW 
This relation is a consequence of four-momentum conservation, 

or equivalently 

p" + = p°' + (6) 

p + q=p' + q'. (7) 

These are the conservation equations for relativistic particle dynamics. In the classical case the 
corresponding conservation equations read 

\p\^ + \q\' = \p'\^ + \q'\\ (8) 
p + q^p' + q'. (9) 

The function a{p, q, uj) gives the distance between p and p' {q and q') in momentum space, and the 
analogue function in the non-relativistic, Newtonian, classical case has the form 

aci(p,g,cj) = w • (q-p). (10) 

By inserting Oci in place of a in Equation ([3| we obtain the classical Boltzmann collision operator 
(disregarding the scattering kernel, which is also different). We point out that there are other 
representations of the colhsion operator ([s]), cf. [179] . 

In [33] and [178) classical solutions to the relativistic Boltzmann equations are studied as c — >■ oo, 
and it is proven that the limit as c — > cx) of these solutions satisfies the classical Boltzmann 
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equation. The former work is more general since general initial data is considered, whereas the 
latter is concerned with data near vacuum. The latter result is stronger in the sense that the limit, 
as c — ^ oo, is shown to be uniform in time. 

The main result concerning the existence of solutions to the classical Boltzmann equation is 
a theorem by DiPerna and Lions [71] that proves existence, but not uniqueness, of renormalized 
solutions. An analogous result holds in the relativistic case, as was shown by Dudyhski and 
Ekiel-Jezewska [75], cf. also |102| . Regarding classical solutions, lUner and Shinbrot [99^ have 
shown global existence of solutions to the non-relativistic Boltzmann equation for initial data close 
to vacuum. Glassey showed global existence for data near vacuum in the relativistic case in a 
technical work |80j . He only requires decay and integrability conditions on the differential cross- 
section, although these are not fully satisfactory from a physics point of view. By imposing more 
restrictive cut-off assumptions on the differential cross-section. Strain [178 gives a different proof, 
which is more related to the proof in the non-relativistic case [99j than [8C), is. For the homogeneous 
relativistic Boltzmann equation, global existence for small initial data has been shown in |12()j under 
the assumption of a bounded differential cross-section. For initial data close to equilibrium, global 
existence of classical solutions has been proven by Glassey and Strauss [87j using assumptions on 
the differential cross-section, which fall into the regime "hard potentials" , whereas Strain |177j has 
shown existence in the case of soft potentials. Rates of the convergence to equilibrium are given in 
both [57] and |177j . In the non-relativistic case, we refer to [1891 1172[ 1119] for analogous results. 

The collision operator Q{f,g) may be written in an obvious way as 

Q{f,9) = Q+{f,9)-Q-{f,9), 

where (5+ and are called the gain and loss term, respectively. If the loss term is deleted the 
gain-term-only Boltzmann equation is obtained. It is interesting to note that the methods of proof 
for the small data results mentioned above concentrate on gain-term-only equations, and once that 
is solved it is easy to include the loss term. In ^1] it is shown that the gain-term-only classical and 
relativistic Boltzmann equations blow up for initial data not restricted to a small neighborhood 
of trivial data. Thus, if a global existence proof of classical solutions for unrestricted data will be 
given, it will necessarily use the full collision operator. 

The gain term has a nice regularizing property in the momentum variable. In ^ it is proven 
that given / e L'^(M.^) and g e ^^(IR^) with f,g>0, then 

l!Q^(/,5)ll_f/i(R?) < C'll/IU2(K|)ll5llLi(Kg)' (11) 

under some technical requirements on the scattering kernel. Here iJ* is the usual Sobolev space. 
This regularizing result was first proven by Lions [112] in the classical situation. The proof relies 
on the theory of Fourier integral operators and on the method of stationary phase, and requires a 
careful analysis of the collision geometry, which is very different in the relativistic case. Simplified 
proofs in the classical and relativistic case are given in [1931 1194j . 

The regularizing theorem has many applications. An important application is to prove that 
solutions tend to equilibrium for large times. More precisely, Lions used the regularizing theorem 
to prove that solutions to the classical Boltzmann equation, with periodic boundary conditions, 
converge in to a global Maxwellian, 

M = e-"IPl'+'3-P+T with a,-feR, a > 0, /3 € M^ 

as time goes to infinity. This result was first obtained by Arkcryd [29 by using non-standard 
analysis. It should be pointed out that the convergence takes place through a sequence of times 
tending to infinity and it is not known whether the limit is unique or depends on the sequence. In 
the relativistic situation, the analogous question of convergence to a relativistic Maxwellian, or a 
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Jiittner equilibrium solution, 



J ^ g-aVi+|pP+/3-P+7^ q;, /3, and 7 as above, with a > 

was studied by Glassey and Strauss [53 ISH]- In the periodic case, they proved convergence in 
a variety of function spaces for initial data close to a Jiittner solution. Having obtained the 
regularizing theorem for the relativistic gain term, it is a straightforward task to follow the method 
of Lions and prove convergence to a global Jiittner solution for arbitrary initial data (satisfying 
the natural bounds of finite energy and entropy), which are periodic in the space variables, cf. [4]. 
We also mention that in the non-rclativistic case Desvillettes and Villani [BS] have studied the 
convergence rate to equilibrium in detail. A similar study in the relativistic case has not yet been 
achieved. 

For more information on the relativistic Boltzmann equation on Minkowski space we refer 
to O [Ml HHH ES] and in the non-relativistic case we refer to [T90| . [79 l l53 ] . 



1.2 The Vlasov— Maxwell and Vlasov— Poisson systems 

Let us consider a collision-less plasma, which is a collection of particles for which collisions are 
relatively rare and the interaction is through their charges. For simplicity we assume that the 
plasma consists of one type of particle, although the results below hold for plasmas with several 
particle species. The particle rest mass and the particle charge are normalized to one. In the 
kinetic framework, the most general set of equations for modeling a collision-less plasma is the 
relativistic Vlasov-Maxwell system: 

dtf + v- V,/ + {Eit, x)+vx B{t, x)) • V„/ = (12) 

dtE + j^c\/xB, \/-E = p, (13) 

dtB = -cV X E, V -B = 0. (14) 

The notation follows the one already introduced with the exception that the momenta are now 
denoted by v instead of p. This has become a standard notation in this field. E and B are the 
electric and magnetic fields, and v is the relativistic velocity, 

(15) 



v/TThW 

where c is the speed of light. The charge density p and current j are given by 



fdv, j ^ vfdv. (16) 



Equation ( 12 ) is the relativistic Vlasov equation and Equations (13 14 ) are the Maxwell equations. 

A special case in three dimensions is obtained by considering spherically-symmetric initial data. 
For such data it can be shown that the solution will also be spherically symmetric, and that the 
magnetic field has to be constant. The Maxwell equation V x E ~ —dtB then implies that 
the electric field is the gradient of a potential (p. Hence, in the spherically-symmetric case the 
relativistic Vlasov-Maxwell system takes the form 

dtf + v-V^f + PE{t,x)-V^f = Q, (17) 
E = V0, = p. (18) 

Here P — \, and the constant magnetic field has been set to zero, since a constant field has no 
significance in this discussion. This system makes sense for any initial data, without symmetry 
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constraints, and is called the relativistic Vlasov-Poisson system. Another special case of interest 
is the classical limit, obtained by letting c — cx) in Equations ( [T2| [l3| [m] ) , yielding: 

dtf + v-V,f + PEit,x)-V^f^O, (19) 
E ^ V(/), Acj) = p, (20) 

where /3 = 1. We refer to Schaeffer |166j for a rigorous derivation of this result. This is the Vlasov- 
Poisson system, and /3 = 1 corresponds to repulsive forces (the plasma case). Taking /3 = — 1 means 
attractive forces and the Vlasov-Poisson system is then a model for a Newtonian self-gravitating 
system. 

One of the fundamental problems in kinetic theory is to find out whether or not spontaneous 
shock formations will develop in a collision- less gas, i.e., whether solutions to any of the equations 
above will remain smooth for all time, given smooth initial data. 

If the initial data are small this problem has an affirmative solution in all cases considered 
above [511 IHSl [3H 133] • For initial data unrestricted in size the picture is more involved. In order to 
obtain smooth solutions globally in time, the main issue is to control the support of the momenta 

Q{t) sup{|w| : 3(s, x) e [0, t] x such that /(s, x, v) ^ 0}, (21) 

i.e., to bound Q{t) by a continuous function so that Q{t) will not blow up in finite time. That such 
a control is sufficient for obtaining global existence of smooth solutions follows from well-known 
results in the different cases, cf. [85 l 11041 [39l[96ll34ll8T| . For the full three-dimensional relativistic 
Vlasov-Maxwell system, the problem of establishing whether or not solutions will remain smooth 
for all time is open. A different sufficient criterion for global existence in this case is given by 
Pallard in |129j . and he also shows a new bound for the electromagnetic field in terms of Q{t) 
in '130 . In two space and three momentum dimensions, Glassey and Schaeffer [821 183j have shown 
that Q{t) can be controlled for the relativistic Vlasov-Maxwell system, which thus yields global 
existence of smooth solutions in that case. 

The relativistic and non-relativistic Vlasov-Poisson equations are very similar in form. In 
particular, the equation for the field is identical in the two cases. However, the mathematical 
results concerning the two systems are very different. In the non-relativistic case, Batt [31] gave 
an affirmative solution in 1977 in the case of spherically-symmetric data. Pfaffelmoser |133) was 
the first one to give a proof for general smooth data. A simplified version of the proof is given by 
Schaeffer in |168j . Pfaffelmoser obtained the bound 

Q{t) < C(l+i)(5i+*)/", 

where ^ > can be taken as arbitrarily small. This bound was later improved by different authors. 
The sharpest bound valid for (3—1 and (3 — ~1 has been given by Horst [07] and reads 

Q{t) < C(l + t)log(2 + t). 

In the case of repulsive forces (/? = 1) Rein (l37j has found a better estimate by using a new 
identity for the Vlasov-Poisson system, discovered independently by lUner and Rein [55] and by 
Perthame [132) . Rein's estimate reads 

Q{t) < C{l + tf/^ 

Independently, and at about the same time as Pfaffelmoser gave his proof, Lions and Perthame [113] 
used a different method for proving global existence. Their method is more generally applicable, 
and the two studies and |105j are examples of problems in related systems, where their method 
has been successful. On the other hand, their method does not give such strong growth estimates 
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on Q{t) as described above. For the relativistic Vlasov-Poisson equation, Glassey and Schaeffer [5T] 
showed in the case /3 = 1 that if the data are spherically symmetric, Q{t) can be controlled, which 
is analogous to the result by Batt mentioned above. Also in the case of cylindrical symmetry 
they are able to control Q{t); see [Mj. If /? = — 1 it was shown in that blow-up occurs in 
finite time for spherically-symmetric data with negative total energy. More recently, Lemou et 
al. [Ill] have investigated the structure of the blow-up solution. They show that the blow-up is 
determined by the self-similar solution of the ultra-relativistic gravitational Vlasov-Poisson system. 
It should be pointed out that the relativistic Vlasov-Poisson system is unphysical since it lacks 
the Lorentz invariance; it is a hybrid of a classical Galilei invariant field equation and a relativistic 
transport equation (17), cf. [3]. In particular, in the case /3 = — 1, it is not a special case of the 
Einstein- Vlasov system. Only for spherically-symmetric data, in the case (3 — 1, is the equation 
a fundamental physical equation. The results mentioned above all concern classical solutions. 
The situation for weak solutions is different, in particular the existence of weak solutions to the 
relativistic Vlasov-Maxwell system is known [701 1139j . 

We also mention that models, which take into account both collisions and the electric and mag- 
netic fields generated by the particles have been investigated. Classical solutions near a Maxwellian 
for the Vlasov-Maxwell-Boltzmann system are constructed by Guo in [SU]. A similar result for 
the Vlasov-Maxwell-Landau system near a Jiittner solution is shown by Guo and Strain in |180) . 

We refer to the book by Glassey [73] and the review article by Rein [141] for more information 
on the relativistic Vlasov-Maxwell system and the Vlasov-Poisson system. 



1.3 The Nordstrom— Vlasov system 



Before turning to the main theme of this review, i.e., the Einstein- Vlasov system, we briefly review 
the results on the Nordstrom- Vlasov system. Nordstrom gravity |120j is an alternative theory of 
gravity introduced in 1913. By coupling this model to a kinetic description of matter the Nord- 
strom- Vlasov system results. In Nordstrom gravity the scalar field (f> describes the gravitational 
field in the sense given below. The Nordstrom- Vlasov system reads 



A, 



„40 



fdp 



(22) 
(23) 



Here 



P - 



denotes the relativistic velocity of a particle with momentum p. The mass of each particle, the 
gravitational constant, and the speed of light are all normalized to one. A solution (f, (/>) of this 
system is interpreted as follows. The spacetime is a Lorentzian manifold with a conformally-flat 
metric 

diag(-l, 1,1,1). 



The particle distribution / defined on the mass shell in this metric is given by 

f{t,x,p)^f{t,x,e'^p). 



(24) 



The first mathematical study of this system was initiated by Calogero in [13] , where the existence 
of static solutions is established. The stability of the static solutions was then investigated in [52] . 
Although the Nordstrom- Vlasov model of gravity does not describe physics correctly, the system 
approaches the Vlasov-Poisson system in the classical limit. Indeed, it is shown in that 
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solutions of the Nordstrom- Vlasov system tend to solutions of the Vlasov-Poisson system as the 
speed of light goes to infinity. 

The Cauchy problem was studied by several authors [STJ jSOl [HI 11081 1131| and the question of 
global existence of classical solutions for general initial data was open for some time. The problem 
was given an affirmative solution in 2006 by Calogero |45^. Another interesting result for the 
Nordstrom- Vlasov system is given in |36j , where a radiation formula, similar to the dipole formula 
in electrodynamics, is rigorously derived. 
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2 The Einstein— Vlasov System 



In this section wc consider a self-gravitating collision-less gas in the framework of general relativity 
and we present the Einstein- Vlasov system. It is most often the case in the mathematics literature 
that the speed of light c and the gravitational constant G are normalized to one, but we keep these 
constants in the formulas in this section since in some problems they do play an important role. 
However, in most of the problems discussed in the forthcoming sections these constants will be 
normalized to one. 

Let M be a four-dimensional manifold and let gab be a metric with Lorentz signature (— , 
so that (A/, gab) is a spacetime. The metric is assumed to be time-orientable so that there is a 
distinction between future and past directed vectors. 

The possible values of the four-momentum p° of a particle with rest mass m belong to the mass 
shell Pm C TM, defined by 

P„, ■■= {{x"-,?") G TM : gab{x°')p°'P^ = -m^c^, p" is future directed}. (25) 

Hence, if m > 0, Pmix"") is the set of all future-directed time- like vectors with length cm, and if 
m = it is the set of all future-directed null vectors. On Pm we take {x°',p'), a = 0,1,2,3 and 
i = 1,2,3 (letters in the beginning of the alphabet always take values 0, 1, 2, 3 and letters in the 
middle take 1, 2, 3) as local coordinates, and is expressed in terms of pp and the metric in view 



of Equation (25 1 . Thus, the density function is a non-negative function on P^. Below we drop 
the index m on /,„ and simply write /. 

Since we are considering a collisionless gas, the particles follow the geodesies in spacetime. The 
geodesies are projections onto spacetime of the curves in Pm defined in local coordinates by 



dpi 
ds 



P 



= -ncpy 



Here F^^ are the Christoffel symbols. Along a geodesic the density function / = f{x'^,p^) is 
invariant so that 

ds 

which implies that 

This is accordingly the Vlasov equation. We point out that sometimes the density function is 
considered as a function on the entire tangent bundle TM rather than on the mass shell Pm C TM. 
The Vlasov equation for / = /(x",p°) then takes the form 



This equation follows from (26 1 if we take the mass shell condition gabP°'P — —ra into account 



Indeed, by abuse of notation, we have 



9/ 


^5/ 




dx"^ 


df 


^ a/ 


dp^ 


dpi 
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Here / is considered as a function on P„j in the left-hand side, and on TAI in the right-hand side. 
From the mass shell condition gabP'^p'' = —m?c^ we derive 



1 

Po' 

.El 
Po' 



-p'PcK, 



Inserting these relations into (26) we obtain (27). If we let t — x'^, and divide the Vlasov equation 
( 26 ) by p'^ we obtain the most common form in the literature of the Vlasov equation 



dt pO dx^ pO ^ dp3 



(28) 



In a fixed spacetime the Vlasov equation ( 28 ) is a linear hyperbolic equation for / and we can 
solve it by solving the characteristic system, 



dP' _ : P^P^ 
~di ~~ '^''~p5~- 

In terms of initial data /o the solution of the Vlasov equation can be written as 

f{x%p') = MX\0,x-,f),P\0,x-,p')), 



(29) 
(30) 

(31) 



where X*(s,x'',p*) and P^{s, x°- ,p^) solve Equations (29 30), and where 

X*(^,a;^p*) = x' and P\t,x'',p') = p\ 

In order to write down the Einstein- Vlasov system we need to know the energy-momentum 
tensor in terms of / and Qab- We define 



Vl.9a6| / fPa 



dp^dp^dp^ 



aPb- 



-Po 



(32) 



where, as usual, pa — gabP^, and \gab\ denotes the absolute value of the determinant of gab- We 
remark that the measure 

V\9ab\ 



-dp^dp^dp^ 



-Po 



^ab — 



is the induced metric of the submanifold Pm{x'^) C T^aM, and that /i is invariant under Lorentz 
transformations of the tangent space, and it is often the case in the literature that is written 
as 

/ fPaPbfJ- 

Let us now consider a coUisionless gas consisting of particles with different rest masses 
mi, 1712, . . . , rriN, described by N density functions fm^ , j = 1, . . . , N . Then the Vlasov equations 
for the different density functions fmj , together with the Einstein equations, 

Rab - ^Rgab + Agab = J2 ^ab" ' 

fc=l 
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form the Einstein- Vlasov system for the coUision-less gas. Here Rab is the Ricci tensor, R is the 
scalar curvature and A is the cosmological constant. 

Henceforth, we always assume that there is only one species of particles in the gas and we write 
Tab for its energy momentum tensor. Moreover, in what follows, we normalize the rest mass m 
of the particles, the speed of light c, and the gravitational constant G, to one, if not otherwise 
explicitly stated that this is not the case. 

Let us now investigate the features of the energy momentum tensor for Vlasov matter. We 
define the particle current density 



Using normal coordinates based at a given point and assuming that / is compactly supported, it 
is not hard to see that Tab is divergence-free, which is a necessary compatibility condition since 
the left-hand side of ^ is divergence-free by the Bianchi identities. A computation in normal 
coordinates also shows that is divergence-free, which expresses the fact that the number of 
particles is conserved. The definitions of Tab and N"^ immediately give us a number of inequalities. 
If V"' is a future-directed time- like or null vector then we have NaV^ < with equality if and only 
if / = at the given point. Hence, N"' is always future-directed time-like, if there are particles 
at that point. Moreover, if and are future-directed time-like vectors then TabV°'W^ > 0, 
which is the dominant energy condition. This also implies that the weak energy condition holds. 
If AT" is a space- like vector, then TabX°-X^ > 0. This is called the non-negative pressure condition, 
and it implies that the strong energy condition holds as well. That the energy conditions hold 
for Vlasov matter is one reason that the Vlasov equation defines a well-behaved matter model 
in general relativity. Another reason is the well-posedness theorem by Choquet-Bruhat [SFj for 
the Einstein- Vlasov system that we state below. Before stating that theorem we first discuss the 
conditions imposed on the initial data. 

The initial data in the Cauchy problem for the Einstein- Vlasov system consist of a 3-dimcnsional 
manifold S, a Riemannian metric gij on S, a symmetric tensor kij on S, and a non-negative scalar 
function /o on the tangent bundle TS of S. 

The relationship between a given initial data set (gij , kij ) on S and the metric gab on the 
spacetime manifold is, that there exists an embedding i/) of 5 into the spacetime such that the 
induced metric and second fundamental form of ip{S) coincide with the result of transporting 
(gij, kij) with ■(/;. For the relation of the distribution functions / and /o we have to note that / is 
defined on the mass shell. The initial condition imposed is that the restriction of / to the part of 
the mass shell over iJj{S) should be equal to fo o (■0"^, d{tp)~^) o 0, where 4> sends each point of the 
mass shell over ipiS) to its orthogonal projection onto the tangent space to ^/'(S'). An initial data 
set for the Einstein- Vlasov system must satisfy the constraint equations, which read 



Here p — Tabn°"n^ and — —h°'^Tbc'rf , where is the future directed unit normal vector to the 
initial hypersurface, and h""^ = g""^ -\- n°"n!' is the orthogonal projection onto the tangent space to 
the initial hypersurface. In terms of /o we can express p and by [j"" satisfies Uaj"" = 0, so it can 
naturally be identified with a vector intrinsic to S) 




R-kijk'^ + (tr/c)2 = Wnp, 
\7ikl~\/i{trk)^8TTji. 



(33) 
(34) 
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We can now state the local existence theorem by Choquet-Bruhat [SS], for the Einstein- Vlasov 
system. 



Theorem 1 Let S be a 3-dimensional manifold, gij a smooth Riemannian metrie on S, kij a 
smooth symmetric tensor on S and /q a smooth non-negative function of compact support on the 



tangent bundle TS of S . Suppose that these objects satisfy the constraint equations (33 34-)- Then 
there exists a smooth spacetime {M,gab), a smooth distribution function f on the mass shell of 
this spacetime, and a smooth embedding ip of S into M , which induces the given initial data on S 
such that gab and f satisfy the Einstein-Vlasov system and ip{S) is a Cauchy surface. Moreover, 
given any other spacetime {M' ,g'^^), distribution function f and embedding ^' satisfying these 
conditions, there exists a diffeomorphism x from an open neighborhood of il){S) in M to an open 
neighborhood of ^'{S) in M' , which satisfies x ° "0 = V"' '^^d carries gab cind f to and /', 
respectively. 



The above formulation is in the case of smooth initial data; for information on the regularity 
needed on the initial data wc refer to [5S] and |118j . In this context we also mention that local 
existence has been proven for the Yang-Mills-Vlasov system in , and that this problem for the 
Einstein-Maxwell-Boltzmann system is treated in [50] ■ However, this result is not complete, as the 
non-negativity of / is left unanswered. Also, the hypotheses on the scattering kernel in this work 
leave some room for further investigation. The local existence problem for physically reasonable 
assumptions on the scattering kernel does not seem well understood in the context of the Einstein- 
Boltzmann system, and a careful study of this problem would be desirable. The mathematical study 
of the Einstein-Boltzmann system has been very sparse in the last few decades, although there has 
been some activity in recent years. Since most questions on the global properties are completely 
open let us only very briefly mention some of these works. Mucha jli7j has improved the regularity 
assumptions on the initial data assumed in [3(]]. Global existence for the homogeneous Einstein- 
Boltzmann system in Robertson- Walker spacetimes is proven in |125j . and a generalization to 
Bianchi type I symmetry is established in |124| . 

In the following sections we present results on the global properties of solutions of the Einstein- 
Vlasov system, which have been obtained during the last two decades. 

Before ending this section we mention a few other sources for more background on the Einstein- 
Vlasov system, cf. [TSm 1751 [T7B] . 



14 



3 The Asymptotically-Flat Cauchy Problem: Spherically- 
Symmetric Solutions 

In this section, we discuss results on global existence and on the asymptotic structure of solutions 
of the Cauchy problem in the asymptotically-flat case. 

In general relativity two classes of initial data are distinguished in the study of the Cauchy 
problem: asymptotically-flat initial data and cosmological initial data. The former type of data 
describes an isolated body. The initial hypersurface is topologically E'^ and appropriate fall-off 
conditions are imposed to ensure that far away from the body spacetime is approximately flat. 
Spacetimes, which possess a compact Cauchy hypersurface, are called cosmological spacetimes, 
and data are accordingly given on a compact 3-manifold. In this case, the whole universe is 
modeled rather than an isolated body. 

The symmetry classes that admit asymptotic flatness are few. The important ones are spher- 
ically symmetric and axially symmetric spacetimes. One can also consider a case, which is un- 
physical, in which spacetime is asymptotically flat except in one direction, namely cylindrically- 
symmetric spacetimes, cf. [75], where the Cauchy problem is studied. The majority of the work so 
far has been devoted to the spherically-symmetric case but recently a result on static axisymmetric 
solutions has been obtained. 

In contrast to the asymptotically-flat case, cosmological spacetimes admit a large number of 
symmetry classes. This provides the possibility to study many special cases for which the difficulties 
of the full Einstein equations are reduced. The Cauchy problem in the cosmological case is reviewed 
in Section m 

The following subsections concern studies of the spherically-symmetric Einstein- Vlasov system. 
The main goal of these studies is to provide an answer to the weak and strong cosmic censorship 
conjectures, cf. |191l [CT] for formulations of the conjectures. 



3.1 Set up and choice of coordinates 

The study of the global properties of solutions to the spherically-symmetric Einstein- Vlasov system 
was initiated two decades ago by Rein and Rendall |142j , cf. also jl351 1156j . They chose to work 
in coordinates where the metric takes the form 

where t G M., r > 0, 6 E [0, vr], (p E [0, 2tt]. These are called Schwarzschild coordinates. Asymptotic 
flatness is expressed by the boundary conditions 

lim X{t,r) = lim fj.{t,r) = 0, Vi > 0. 

r— >C30 r— ^oo 

A regular center is also required and is guaranteed by the boundary condition 

A(i,0) = yt>o. 

The coordinates (r, 6, <j>) give rise to difficulties at r = and it is advantageous to use Cartesian 
coordinates. With 

X = {r sin (/) cos 9, r sin (f> sin 9, r cos (f)) 



as spatial coordinates and 



r r 
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as momentum coordinates, the Einstein- Vlasov system reads 

dtf + e^-^^J= ■ VJ - (a, — + e^-^^lrV^+W) - ' V„/ = 0, (35) 
V1 + H2 \ r J r 

e-^^{2rXr -!) + ! = STrr^p, (36) 

e"^^(2r/ir + 1) - 1 = 87rrV (37) 

At = -47rre^+'"j, (38) 

(/irr + (Mr - K)(pr + ")) - e'"^" (Xu + At(Af - fit)) = SlTpT . (39) 



-2A 



r 



The matter quantities are defined by 

p{t,x)= I ^l + \v\^f{t,x,v)dv, (40) 



P(t,x)^f (^Yfit^x^v) (41) 
j(t,a;) = / —f{t,x,v)dv, (42) 



PT{t,x) = l [ \^^\''f{t,x,v)dv. (43) 
Here p is the energy density, j the current, p the radial pressure, and px the tangential pressure. 



Let us point out that these equations are not independent, e.g.. Equations (38) and (39) follow 



from (35)-(37) 



As initial data we take a spherically-symmetric, non-negative, and continuously differentiable 
function /o with compact support that satisfies 

^/l+Wfoiy,v)dvdy (44) 

This condition guarantees that no trapped surfaces are present initially. 

The set up described above is one of several possibilities. The Schwarzschild coordinates have 
the advantage that the resulting system of equations can be written in a quite condensed form. 
Moreover, for most initial data, solutions are expected to exist globally in Schwarzschild time, which 
sometimes is called the polar time gauge. Let us point out here that there are initial data leading to 
spacetime singularities, cf. |149[ [20 l I24j. Hence, the question of global existence for general initial 
data is only relevant if the time slicing of the spacetime is expected to be singularity avoiding, 
which is the case for Schwarzschild time. We refer to [116] for a general discussion on this issue. 
This makes Schwarzschild coordinates tractable in the study of the Cauchy problem. However, 
one disadvantage is that these coordinates only cover a relatively small part of the spacetime, in 
particular trapped surfaces are not admitted. Hence, to analyze the black-hole region of a solution 
these coordinates are not appropriate. Here we only mention the other coordinates and time gauges 
that have been considered in the study of the spherically symmetric Einstein- Vlasov system. These 
works will be discussed in more detail in various sections below. Rendall uses maximal-isotropic 
coordinates in |156j . These coordinates are also considered in [T^]. The Einstein- Vlasov system 
is investigated in double null coordinates in [64l [63] . Maximal-areal coordinates and Eddington- 
Finkelstein coordinates are used in [3TJ |T7j , and in [21] respectively. 

3.2 Local existence and the continuation criterion 



In [142] it is shown that for initial data satisfying (44) there exists a unique, continuously 



differentiable solution / with /(O) = /o on some right maximal interval [0,T). If the solution 
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blows up in finite time, i.e., if T < oo, then p{t) becomes unbounded as t — > T. Moreover, a 
continuation criterion is shown that says that a local solution can be extended to a global one, 
provided Q{t) can be bounded on [0,r), where 

Q{t) sup{|w| : 3(s, x) G [0, t] x such that /(s, x, v) ^ 0}. (45) 

This is analogous to the situation for the Vlasov-Maxwell system. A control of the w-support 



immediately implies that p and p are bounded in view of Equations (40 41 ). In the Vlasov-Maxwell 
case the field equations have a regularizing effect in the sense that derivatives can be expressed 
through spatial integrals, and it follows [85] that the derivatives of / can also be bounded if the 
ii-support is bounded. For the Einstein- Vlasov system such a regularization is less clear, since, e.g.. 



Iir depends on p in a point-wise manner. However, in view of Equation (39) certain combinations 
of second and first order derivatives of the metric components can be expressed in terms of the 
matter component pT, which is a consequence of the geodesic deviation equation. This fact turns 
out to also be sufficient for obtaining bounds on the derivatives of /, cf. |1421I135T I156j for details. 

The local existence result discussed above holds for compactly-supported initial data. The 
compact support condition in the momentum variables is in ^2, replaced by the fall-off condition 

sup {l + \v\f\fix,v)\<cx,. (46) 

We also refer to [53] where a subclass of non-compactly-supported data is treated. 

Local existence of solutions in double null coordinates and in Eddington-Finkelstein coordinates 
is established in [B^, and [2^ respectively. 



3.3 Global existence for small initial data 

In [142j the authors also consider the problem of global existence in Schwarzschild coordinates for 
small initial data for massive particles. They show that for such data the w-support is bounded on 
[0,T). Hence, the continuation criterion implies that T = oo. The resulting spacetime in [142] is 
geodesically complete, and the components of the energy-momentum tensor as well as the metric 
quantities decay with certain algebraic rates in t. The mathematical method used by Rein and 
Rendall is inspired by the analogous small data result for the Vlasov-Poisson equation by Bardos 
and Degond ^32j. This should not be too surprising since for small data the gravitational fields are 
expected to be small and a Newtonian spacetime should be a fair approximation. In this context 
we point out that in [143] it is proven that the Vlasov-Poisson system is indeed the non-relativistic 
limit of the spherically-symmetric Einstein- Vlasov system, i.e., the limit when the speed of light 
c — >■ cxD. In |150j this result is shown without symmetry assumptions. 

As mentioned above the local and global existence problem has been studied using other time 
gauges, in particular Rendall has shown global existence for small initial data in maximal-isotropic 
coordinates in |156j . 

The previous results refer to massive particles but they do not immediately carry over to 
massless particles. This case is treated by Dafermos in [63] where global existence for small initial 
data is shown in double null coordinates. The spacetimes obtained in the studies [1421 1156( [63] 
are all causally geodesically complete and appropriate decay rates of the metric and the matter 
quantities are given. 



3.4 Global existence for special classes of large initial data 

In the case of small initial data the resulting spacetime is geodesically complete and no singularities 
form. A different scenario, which leads to a future geodesically complete spacetime, is to consider 
initial data where the particles are moving rapidly outwards. If the particles move sufficiently fast 
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the matter disperses and the gravitational attraction is not strong enough to reverse the velocities 
of the particles to create a collapsing system. This problem is studied in [T7j using a maximal 
time coordinate. It is shown that the scenario described above can be realized, and that global 
existence holds. 

In Section |3.7| we discuss results on the formation of black holes and trapped surfaces; in 
particular, the results in [20] will be presented. A corollary of the main result in [20 concerns 
the issue of global existence and thus we mention it here. It is shown that a particular class of 
initial data, which lead to formation of black holes, have the property that the solutions exist 
for all Schwarzschild time. The initial data consist of two parts: an inner part, which is a static 
solution of the Einstein- Vlasov system, and an outer part with matter moving inwards. The set-up 
is shown to preserve the direction of the momenta of the outer part of the matter, and it is also 
shown that in Schwarzschild time the inner part and the outer part of the matter never interact 
in Schwarzschild time. 



3.5 On global existence for general initial data 



As was mentioned at the end of Section 3.1 the issue of global existence for general initial data 
is only relevant in certain time gauges since there are initial data leading to singular spacetimes. 
However, it is reasonable to believe that global existence for general data may hold in a polar time 
gauge or a maximal time gauge, cf. 1116] . and it is often conjectured in the literature that these time 
slicings are singularity avoiding. However, there is no proof of this statement for any matter model 
and it would be very satisfying to provide an answer to this conjecture for the Einstein- Vlasov 
system. A proof of global existence in these time coordinates would also be of great importance 
due to its relation to the weak cosmic censorship conjecture, cf. [61 ] l62 | l65]. 

The methods of proofs in the cases described in Sections |3.3| and |3.4[ where global existence 
has been shown, are all tailored to treat special classes of initial data and they will likely not apply 
in more general situations. In this section we discuss some attempts to treat general initial data. 
These results are all conditional in the sense that assumptions are made on the solutions, and not 
only on the initial data. 

The first study on global existence for general initial data is [146j . which is carried out in 
Schwarzschild coordinates. The authors introduce the following variables in the momentum space 
adapted to spherical symmetry, 

L:^\x\''\v\^ -{x-vf, w^—, (47) 

r 

where L is the square of the angular momentum and w is the radial component of the momenta. 
A consequence of spherical symmetry is that angular momentum is conserved along the character- 
istics. In these variables the Vlasov equation for / = f{t, r, w, L) becomes 

dtf + e^'-^'^drf - (^Xtw + e^^-^iirE - e^"^^) d^f = 0, (48) 

where 

E = E{r, w, L) = yjl + w"^ + L/r"^. 
The main result in |146j shows that as long as there is no matter in the ball 

{x G : |a;| < e}, 

the estimate 

Q(i) <eiogQ(o)eC(^)'^ (49) 
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holds. Here C(e) is a constant, which depends on e. Thus, in view of the continuation criterion 
this can be viewed as a global existence result outside the center of symmetry for initial data with 
compact support. This result rules out shell-crossing singularities, which are present when, e.g., 
dust is used as a matter model. The bound of Q is obtained by estimating each term individually 



in the characteristic equation associated with the Vlasov equation ( 48 ) for the radial momentum. 



This involves a particular difficulty. The Einstein equations imply that 



TO 



where 



m(t, r) = 47r / rj'^ p{t,ri) dr], (50) 







is the quasi local mass. Thus, using (381 the characteristic equation consists of the two terms 
Ti = 4:77 re'^^-^{jw + pE), and T2 — e'^^^^^, together with a term, which is independent of the 
matter quantities. There is a distinct difference between the terms Ti and T2 due to the fact that 
TO can be regarded as an average, since it is given as a space integral of the energy density p, 
whereas j and p are point-wise terms. The method in |146| makes use of a cancellation property 
of the radial momenta in Ti so that outside the center this term is manageable but in general it 
seems very unpleasant to have to treat point-wise terms of this kind. 

In |156| Rendall shows global existence outside the center in maximal-isotropic coordinates. 
The bound on Q{t) is again obtained by estimating each term in the characteristic equation. In 
this case there are no point-wise terms in contrast to the case with Schwarzschild coordinates. 
However, the terms are, in analogy with the Schwarzschild case, strongly singular at the center. 

A recent work [12 gives an alternative and simplified proof of the result in [146j. In particular, 
the method avoids the point-wise terms by using the fact that the characteristic system can be 
written in a form such that Green's formula in the plane can be applied. This results in a com- 
bination of terms involving second-order derivatives, which can be substituted for by one of the 
Einstein equations. This method was first introduced in T] but the set-up is different in and 
the application of Green's formula becomes very natural. In addition, the bound of Q is improved 



compared to ( 49 ) and reads 

Q{t) < g(0)e^(i+*)/^'. 

This bound is sufficient to conclude that global existence outside the center also holds for non- 
compact initial data. In addition to the global existence result outside the centre, it is shown in [T^] 
that as long as 3m{t^ r) < r and j < 0, singularities cannot form. Note that in Schwarzschild coor- 
dinates 2m{t, r) < r always, and that there are closed null geodesies if 3to = r in the Schwarzschild 
static spacetime. 

The method in [T^ also applies to the case of maximal-isotropic coordinates studied in |156) . 
There is an improvement concerning the regularity of the terms that need to be estimated to obtain 
global existence in the general case. A consequence of [T^] is accordingly that the quite different 
proofs in |146' and in |156j are put on the same footing. We point out that the method can also 
be applied to the case of maximal-areal coordinates. 

The results discussed above concern time gauges, which are expected to be singularity avoiding 
so that the issue of global existence makes sense. An interpretation of these results is that "first 
singularities" (where the notion of "first" is tied to the causal structure), in the non-trapped 
region, must emanate from the center and that this case has also been shown in double null- 
coordinates by Dafermos and Rendall in [64j . The main motivation for studying the system in 
these coordinates has its origin from the method of proof of the cosmic-censorship conjecture for 
the Einstein-scalar field system by Christodoulou [30] • An essential part of his method is based 
on the understanding of the formation of trapped surfaces [35]. In '55] it is shown that a single 
trapped surface or marginally-trapped surface in the maximal development implies that weak 
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cosmic censorship holds The theorem holds true for any spherically-symmetric matter spacetime if 
the matter model is such that "first" singularities necessarily emanate from the center. The results 
in |146j and in |156j are not sufficient for concluding that the hypothesis of the matter needed in 
the theorem in |62j is satisfied, since they concern a portion of the maximal development covered 
by particular coordinates. Therefore, Dafermos and Rendall [64] choose double-null coordinates, 
which cover the maximal development, and they show that the mentioned hypothesis is satisfied 
for Vlasov matter. 



3.6 Self-similar solutions 

The main reason that the question of global existence in certain time coordinates discussed in the 



previous Section 3.5 is of great importance is its relation to the cosmic censorship conjectures. 
Now there is, in fact, no theorem in the literature, which guarantees that weak cosmic censorship 
follows from such a global existence result, but there are strong reasons to believe that this is the 
case, cf. ,116j and jlSj. Hence, if initial data can be constructed, which lead to naked singularities, 
then either the conjecture that global existence holds generally is false or the viewpoint that global 
existence implies the absence of naked singularities is wrong. In view of a recent result by Rendall 



and Velazquez [161) on self similar dust-like solutions for the massless Einstein- Vlasov system, this 
issue has much current interest. Let us mention here that there is a previous study on self-similar 
solutions in the massless case by Marti'n-Garcfa and Gundlach |115| . However, this result is based 
on a scaling of the density function itself and therefore makes the result less related to the Cauchy 
problem. Also, their proof is, in part, based on numerics, which makes it harder to judge the 
relevance of the result. 

The main aim of the work |161) is to establish self-similar solutions of the massive Einstein- 
Vlasov system and the present result can be viewed as a first step to achieving this. In the set-up, 
two simplifications are made. First, the authors study the massless case in order to find a scaling 
group, which leaves the system invariant. More precisely, the massless system is invariant under 
the scaling 

r ^ Or, t ^ Ot, w — > — L OL. 

The massless assumption seems not very restrictive since, if a singularity forms, the momenta 
will be large and therefore the influence of the rest mass of the particles will be negligible, so 
that asymptotically the solution can be self-similar also in the massive case, cf. for the 

relativistic Vlasov-Poisson system. The second simplification is that the possible radial momenta 
are restricted to two values, which means that the density function is a distribution in this variable. 
Thus, the solutions can be thought of as intermediate between smooth solutions of the Einstein- 
Vlasov system and dust. 

For this simplified system it turns out that the existence question of self-similar solutions can 
be reduced to that of the existence of a certain type of solution of a four-dimensional system of 
ordinary differential equations depending on two parameters. The proof is based on a shooting 
argument and involves relating the dynamics of solutions of the four-dimensional system to that 
of solutions of certain two- and three-dimensional systems obtained from it by limiting processes. 
The reason that an ODE system is obtained is due to the assumption on the radial momenta, 
and if regular initial data is considered, an ODE system is not sufficient and a system of partial 
differential equations results. 

The self-similar solution obtained by Rendall and Velazquez has some interesting properties. 
The solution is not asymptotically flat but there are ideas outlined in |161| of how this can be over- 
come. It should be pointed out here that a similar problem occurs in the work by Christodoulou [SI] 
for a scalar field, where the naked singularity solutions are obtained by truncating self-similar data. 
The singularity of the self-similar solution by Rendall and Velazquez is real in the sense that the 
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Kretschmann scalar curvature blows up. The asymptotic structure of the solution is striking in 
view of the conditional global existence result in [12 . Indeed, the self similar solution is such that 
j < 0, and 3m{t,r) — ^ r asymptotically, but for any T, 3m{t,r) > r for some t > T. In [T^] global 
existence follows if j < and if 3m{t, r) < r for all t. It is also the case that if m/r is close to 1/2, 
then global existence holds in certain situations, cf. [20]. Hence, the asymptotic structure of the 
self-similar solution has properties, which have been shown to be difficult to treat in the search for 
a proof of global existence. 

3.7 Formation of black holes and trapped surfaces 

We have previously mentioned that there exist initial data for the spherically-symmetric Einstein- 
Vlasov system, which lead to formation of black holes. 

The first result in this direction was obtained by Rendall |149j . He shows that there exist initial 
data for the spherically-symmetric Einstein- Vlasov system such that a trapped surface forms in 
the evolution. The occurrence of a trapped surface signals the formation of an event horizon. As 
mentioned above, Dafermos |62| has proven that, if a spherically-symmetric spacetime contains a 
trapped surface and the matter model satisfies certain hypotheses, then weak cosmic censorship 
holds true. In |64] it was then shown that Vlasov matter does satisfy the required hypotheses. 
Hence, by combining these results it follows that initial data exist, which lead to gravitational 
collapse and for which weak cosmic censorship holds. However, the proof in |149j rests on a 
continuity argument, and it is not possible to tell whether or not a given initial data set will 
give rise to a black hole. Moreover, the mechanism of how trapped surfaces form is not revealed 
in |149j . This is in contrast to the result in ;24j, where explicit conditions on the initial data are 
given, which guarantee the formation of trapped surfaces in the evolution. The analysis is carried 
out in Eddington-Finkelstein coordinates and a central result in [21] is to control the life span of 
the solution to ensure that there is sufficient time to form a trapped surface before the solution 
may break down. In particular, weak cosmic censorship holds for these initial data. In |20| the 
formation of the event horizon in gravitational collapse is analyzed in Schwarzschild coordinates. 
Note that these coordinates do not admit trapped surfaces. The initial data in [20] consist of two 
separate parts of matter. One inner part and one outer part, in which all particles move inward 
initially. The reason for the inner part is that it is possible to choose the parameters for the data 
such that the particles of the outer matter part continue to move inward for all Schwarzschild time 
as long as the particles do not interact with the inner part. This fact simplifies the analysis since 
the dynamics is much restricted when the particles keep the direction of their radial momenta. 
The main result is that explicit conditions on the initial data with ADM mass M are given such 
that there is a family of outgoing null geodesies for which the area radius r along each geodesic is 
bounded by 2M. It is furthermore shown that if 

t > 0, and r > 2M + ae~^\ 

where a and P are positive constants, then /(i, r, •, •) = 0, and the metric equals the Schwarzschild 
metric 

ds"" - -(1 - + (1 - "^y'dr^ + r^ide^ + sinHdcf'), (51) 

representing a black hole with mass M . Hence, spacetime converges asymptotically to the Schwarz- 
schild metric. 

The latter result does not reveal whether or not all the matter crosses r = 2M or simply piles 
up at the event horizon. In |23| it is shown that for initial data, which are closely related to 
those in [5Dj, but such that the radial momenta are unbounded, all the matter do cross the event 
horizon asymptotically in Schwarzschild time. This is in contrast to what happens to freely-falling 
observers in a static Schwarzschild spacetime, since they will never reach the event horizon. 
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The result in [20j is reconsidered in |19j . where an additional argument is given to match the 
definition of weak cosmic censorship given in |61j . 

It is natural to relate the results of [^UJ [53] to those of Christodoulou on the spherically- 
symmetric Einstein-scalar field system [57] and [SS]. In [57j it is shown that if the final Bondi 
mass M is different from zero, the region exterior to the sphere r = 2M tends to the Schwarzschild 
metric with mass M similar to the result in [20]. In [SS| explicit conditions on the initial data are 
specified, which guarantee the formation of trapped surfaces. This paper played a crucial role in 
Christodoulou's proof [BD] of the weak and strong cosmic censorship conjectures. The conditions 
on the initial data in [SS] allow the ratio of the Hawking mass and the area radius to cover the full 
range, i.e., 2m/r G (0,1), whereas the conditions in J24j require 2m/r to be close to one. Hence, 
it would be desirable to improve the conditions on the initial data in [24] , although the conditions 
by Christodoulou for a scalar field are not expected to be sufficient in the case of Vlasov matter. 

3.8 Numerical studies on critical collapse 

In |147) a numerical study on critical collapse for the Einstein- Vlasov system was initiated. A 
numerical scheme originally used for the Vlasov-Poisson system was modified to the spherically- 
symmetric Einstein- Vlasov system. It has been shown by Rein and Rodewis |148j that the numer- 
ical scheme has desirable convergence properties. (In the Vlasov-Poisson case, convergence was 
proven in [167] . see also [77]). 

The speculation discussed above that there may be no naked singularities formed for any regular 
initial data is in part based on the fact that the naked singularities that occur in scalar field collapse 
appear to be associated with the existence of type II critical collapse, while Vlasov matter is of 
type I. The primary goal in [147] was indeed to decide whether Vlasov matter is type I or type II. 

These different types of matter are defined as follows. Given small initial data, no black holes 
form and matter will disperse. For large data, black holes will form and consequently there is a 
transition regime separating dispersion of matter and formation of black holes. If we introduce a 
parameter A on the initial data such that for small A dispersion occurs and for large A a black 
hole is formed, we get a critical value A^ separating these regions. If we take A> Ac and denote 
by me (A) the mass of the black hole, then if tob(^) —J' as A —>■ Ac we have type II matter, 
whereas for type I matter this limit is positive and there is a mass gap. For more information on 
critical collapse we refer to the review paper by Gundlach [55] . 

The conclusion of [147 1 is that Vlasov matter is of type I. There are two other independent nu- 
merical simulations on critical collapse for Vlasov matter [1281121] . In these simulations, maximal- 
areal coordinates are used rather than Schwarzschild coordinates as in [147] . The conclusion of 
these studies agrees with the one in [147] . 

3.9 The charged case 

We end this section with a discussion of the spherically-symmetric Einstein- Vlasov-Maxwell sys- 
tem, i.e., the case considered above with charged particles. Whereas the constraint equations in 
the uncharged case, written in Schwarzschild coordinates, do not involve solving any difficulties 
once the distribution function is given, the charged case is more challenging. However, in [123] it is 
shown that solutions to the constraint equations do exist for the Einstein- Vlasov-Maxwell system. 
In [122] local existence is shown together with a continuation criterion, and global existence for 
small initial data is shown in [121] . 
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4 The Cosmological Cauchy Problem 



In this section we discuss the Einstein- Vlasov system for cosmological spacetimes, i.e., spacetimes 
that possess a compact Cauchy surface. The "particles" in the kinetic description are in this case 
galaxies or even clusters of galaxies. The main goal is to determine the global properties of the 
solutions to the Einstein- Vlasov system for initial data given on a compact 3-manifold. In order 
to do so, a global time coordinate t must be found and the asymptotic behavior of the solutions 
when t tends to its limiting values has to be analyzed. This might correspond to approaching a 
singularity, e.g., the big bang singularity, or to a phase of unending expansion. 

Presently, the aim of most of the studies of the cosmological Cauchy problem has been to show 
existence for unrestricted initial data and the results that have been obtained are in cases with 
symmetry (see, however, |27j . where to some extent global properties are shown in the case without 
symmetry). These studies will be reviewed below. A recent and very extensive work by Ringstrom 
has, on the other hand, a different aim, i.e., to show stability of homogeneous cosmological models, 
and concerns the general case without symmetry. The size of the Cauchy data is in this case very 
restricted but, since Ringstrom allows general perturbations, there are no symmetries available to 
reduce the complexity of the Einstein- Vlasov system. This result will be reviewed at the end of 
this section. 

4.1 Spatially-homogeneous spacetimes 

The only spatially-homogeneous spacetimes admitting a compact Cauchy surface are the Bianchi 
types I, IX and the Kantowski-Sachs model; to allow for cosmological solutions with more general 
symmetry types, it is enough to replace the condition that the spacetime is spatially homogeneous, 
with the condition that the universal covering of spacetime is spatially homogeneous. Spacetimes 
with this property are called locally spatially homogeneous and these include, in addition, the 
Bianchi types II, III, V, VIo, VIIq, and VIII. 

One of the first studies on the Einstein- Vlasov system for spatially-homogeneous spacetimes is 
the work [152] by Rendall. He chooses a Gaussian time coordinate and investigates the maximal 
range of this time coordinate for solutions evolving from homogeneous data. For Bianchi IX and 
for Kantowski-Sachs spacetimes he finds that the range is finite and that there is a curvature 
singularity in both the past and the future time directions. For the other Bianchi types there is a 
curvature singularity in the past, and to the future spacetime is causally geodesically complete. In 
particular, strong cosmic censorship holds in these cases. 

Although the questions on curvature singularities and geodesic completeness are very important, 
it is also desirable to have more detailed information on the asymptotic behavior of the solutions, 
and, in particular, to understand in which situations the choice of matter model is essential for the 
asymptotics. 

In recent years several studies on the Einstein- Vlasov system for spatially locally homogeneous 
spacetimes have been carried out with the goal to obtain a deeper understanding of the asymptotic 
structure of the solutions. Roughly, these investigations can be divided into two cases: (i) studies 
on non-locally rotationally symmetric (non-LRS) Bianchi I models and (ii) studies of LRS Bianchi 
models. 

In case (i) Rendall shows in [153 that solutions converge to dust solutions for late times. Under 
the additional assumption of small initial data this result is extended by Nungesser [127] , who gives 
the rate of convergence of the involved quantities. In |153) Rendall also raises the question of the 
existance of solutions with complicated oscillatory behavior towards the initial singularity may 
exist for Vlasov matter, in contrast to perfect fluid matter. Note that for a perfect fluid the 
pressure is isotropic, whereas for Vlasov matter the pressure may be anisotropic, and this fact 
could be sufficient to drastically change the dynamics. This question is answered in [93] . where 
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the existence of a heteroclinic network is established as a possible asymptotic state. This implies 
a complicated oscillating behavior, which differs from the dynamics of perfect fluid solutions. The 
results in |93j were then put in a more general context by Calogero and Heinzle |46) . where quite 
general anisotropic matter models are considered. 

In case (ii) the asymptotic behaviour of solutions has been analyzed in |159[ I160|. l48l l47] . 
In 159", the case of massless particles is considered, whereas the massive case is studied in |160) . 
Both the nature of the initial singularity and the phase of unlimited expansion are analyzed. The 
main concern in these two works is the behavior of Bianchi models I, II, and III. The authors 
compare their solutions with the solutions to the corresponding perfect fluid models. A general 
conclusion is that the choice of matter model is very important since, for all symmetry classes 
studied, there are differences between the collision-less model and a perfect fluid model, both 
regarding the initial singularity and the expanding phase. The most striking example is for the 
Bianchi II models, where they find persistent oscillatory behavior near the singularity, which is 
quite different from the known behavior of Bianchi type II perfect fluid models. In |160| it is also 
shown that solutions for massive particles are asymptotic to solutions with massless particles near 
the initial singularity. For Bianchi I and II, it is also proven that solutions with massive particles 
are asymptotic to dust solutions at late times. It is conjectured that the same also holds true 
for Bianchi III. This problem is then settled by Rendall in [157] . The investigation [l^ concerns 
a large class of anisotropic matter models, and, in particular, it is shown that solutions of the 
Einstein- Vlasov system with massless particles oscillate in the limit towards the past singularity 
for Bianchi IX models. This result is extended to the massive case in 47]. 

Before finishing this section we mention two other investigations on homogeneous models with 
Vlasov matter. In |106j Lee considers the homogeneous spacetimes with a cosmological constant 
for all Bianchi models except Bianchi type IX. She shows global existence as well as future causal 
geodesic completeness. She also obtains the time decay of the components of the energy momentum 
tensor as f — cx), and she shows that spacetime is asymptotically dust-like. Anguige [28] studies the 
conformal Einstein- Vlasov system for massless particles, which admit an isotropic singularity. He 
shows that the Cauchy problem is well posed with data consisting of the limiting density function 
at the singularity. 

4.2 Inhomogeneous models with symmetry 

In the spatially homogeneous case the metric can be written in a form that is independent of the 
spatial variables and this leads to an enormous simplification. Another class of spacetimes that are 
highly symmetric but require the metric to be spatially dependent are those that admit a group of 
isometries acting on two-dimensional spacelike orbits, at least after passing to a covering manifold. 
The group may be two-dimensional (local U{1) x U{1) or symmetry) or three-dimensional 
(spherical, plane, or hyperbolic symmetry). In all these cases, the quotient of spacetime by the 
symmetry group has the structure of a two-dimensional Lorentzian manifold Q. The orbits of 
the group action (or appropriate quotients in the case of a local symmetry) are called surfaces of 
symmetry. Thus, there is a one-to-one correspondence between surfaces of symmetry and points 
of Q. There is a major difference between the cases where the symmetry group is two- or three- 
dimensional. In the three-dimensional case no gravitational waves are admitted, in contrast to 
the two-dimensional case where the evolution part of the Einstein equations are non-linear wave 
equations. 

Three types of time coordinates that have been studied in the inhomogeneous case are CMC, 
areal, and conformal coordinates. A CMC time coordinate t is one where each hypersurface of 
constant time has constant mean curvature and on each hypersurface of this kind the value of t is 
the mean curvature of that slice. In the case of areal coordinates, the time coordinate is a function 
of the area of the surfaces of symmetry, e.g., proportional to the area or proportional to the square 
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root of the area. In the case of conformal coordinates, the metric on the quotient manifold Q is 
conformaUy flat. The CMC and the areal coordinate foHations are both geometrically-based time 
foliations. The advantage with a CMC approach is that the definition of a CMC hypersurface does 
not depend on any symmetry assumptions and it is possible that CMC foliations will exist for 
general spacetimes. The areal coordinate foliation, on the other hand, is adapted to the symmetry 
of spacetime but it has analytical advantages and detailed information about the asymptotics can 
be derived. The conformal coordinates have mainly served as a useful framework for the analysis 
to obtain geometrically-based time foliations. 

4.2.1 Surface symmetric spacetimes 

Let us now consider spacetimes (M, g) admitting a three-dimensional group of isometrics. The 
topology of M is assumed to be M x S"^ x F, with F a compact two-dimensional manifold. The 
universal covering F of F induces a spacetime {M,g) hy M ~ M. x x F and g — p*g, where 
p : M ^ M is the canonical projection. A three-dimensional group G of isometrics is assumed to 
act on {M,g). If F = 5^ and G = SO{3), then iM,g) is called spherically symmetric; if F = 
and G ^ E2 (Euclidean group), then {M,g) is called plane symmetric; and if F has genus greater 
than one and the connected component of the symmetry group G of the hyperbolic plane acts 
isometrically on _F = H^, then {M,g) is said to have hyperbolic symmetry. 

In the case of spherical symmetry the existence of one compact CMC hypersurface implies that 
the whole spacetime can be covered by a CMC time coordinate that takes all real values |151l . 
The existence of one compact CMC hypersurface in this case was proven by Henkel [94] using the 
concept of prescribed mean curvature (PMC) foliation. Accordingly, this gives a complete pic- 
ture in the spherically symmetric case regarding CMC foliations. In the case of areal coordinates, 
Rein [136] has shown, under a size restriction on the initial data, that the past of an initial hyper- 
surface can be covered, and that the Kretschmann scalar blows up. Hence, the initial singularity 
for the restricted data is both a crushing and a curvature singularity. In the future direction it is 
shown that areal coordinates break down in finite time. 

In the case of plane and hyperbolic symmetry, global existence to the past was shown by 
Rendall |151j in CMC time. This implies that the past singularity is a crushing singularity since 
the mean curvature blows up at the singularity. Also in these cases Rein showed |136j under a size 
restriction on the initial data, that global existence to the past in areal time and blow up of the 
Kretschmann scalar curvature as the singularity is approached. Hence, the singularity is both a 
crushing and a curvature singularity in these cases too. In both of these works the question of global 
existence to the future was left open. This gap was closed in [55] , and global existence to the future 
was established in both CMC and areal time coordinates. The global existence result for CMC time 
is a consequence of the global existence theorem in areal coordinates, together with a theorem by 
Henkel [94] which shows that there exists at least one hypersurface with (negative) constant mean 
curvature. In addition, the past direction is analyzed in |25] using areal coordinates, and global 
existence is shown without a size restriction on the data. It is not concluded if the past singularity, 
without the smallness condition on the data, is a curvature singularity as well. The issues discussed 
above have also been studied in the presence of a cosmological constant, cf. |1821I184] . In particular 
it is shown that in the spherically-symmetric case, if A > 0, global existence to the future holds 
in areal time for some special classes of initial data, which is in contrast to the case with A = 0. 
In this context we also mention that surface symmetric spacetimes with Vlasov matter and with a 
Maxwell field have been investigated in |183j . 

An interesting question, which essentially was left open in the studies mentioned above, is 
whether the areal time coordinate, which is positive by definition, takes all values in the range 
(0,00) or only in (i?o,oo) for some positive Rq. It should here be pointed out that there is an 
example for vacuum spacetimes with symmetry (which includes the plane symmetric case) 
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where Rq > 0. This question was first resolved by Weaver |192j for symmetric spacetimes with 
Vlasov matter. Her result shows that if spacetime contains Vlasov matter (/ 7^ 0) then Rq — 0. 
Smulevici [174j has recently shown, under the same assumption, that Rq = also in the hyperbolic 
case. Smulevici also includes a cosmological constant A and shows that both the results, for plane 
(or T"^) symmetry and hyperbolic symmetry, are valid for A > 0. 

The important question of strong cosmic censorship for surface-symmetric spacetimes has re- 
cently been investigated by neat methods by Dafermos and Rendall [571155] . The standard strategy 
to show cosmic censorship is to either show causal geodesic completeness in case there are no singu- 
larities, or to show that some curvature invariant blows up along any incomplete causal geodesic. 
In both cases no causal geodesic can leave the maximal Cauchy development in any extension 
if we assume that the extension is C^. In [67l |65] two alternative approaches are investigated. 
Both of the methods rely on the symmetries of the spacetime. The first method is independent of 
the matter model and exploits a rigidity property of Cauchy horizons inherited from the Killing 
fields. The areal time described above is defined in terms of the Killing fields and a consequence 
of the method by Dafermos and Rendall is that the Killing fields extend continuously to a Cauchy 
horizon, if one exists. Now, since global existence has been shown in areal time it follows that 
there cannot be an extension of the maximal hyperbolic development to the future. This method 
is useful for the expanding future direction. The second method is dependent on Vlasov matter 
and the idea is to follow the trajectory of a particle, which crosses the Cauchy horizon and shows 
that the conservation laws for the particle motion associated with the symmetries of the spacetime, 
such as the angular momentum, lead to a contradiction. In most of the cases considered in [67J 
there is an assumption on the initial data for the Vlasov equation, which implies that the data 
have non-compact support in the momentum space. It would be desirable to relax this assump- 
tion. The results of the studies [671 165| can be summarized as follows. For plane and hyperbolic 
symmetry strong cosmic censorship is shown when A > 0. The restriction that matter has non- 
compact support in the momentum space is here imposed except in the plane case with A = 0. In 
the spherically-symmetric case cosmic censorship is shown when A = 0. In the case of A > a 
detailed geometric characterization of possible boundary components of spacetime is given. The 
diflaculties to show cosmic censorship in this case are related to possible formation of extremal 
Schwarzschild-de-Sitter-type black holes. Cosmic censorship in the past direction is also shown for 
all symmetry classes, and for all values of A, for a special class of anti-trapped initial data. 

Although the methods developed in [67l |65] provide a lot of information on the asymptotic 
structure of the solutions, questions on geodesic completeness and curvature blow up are not 
answered. In a few cases, information on these issues has been obtained. As mentioned above, 
blow up of the Kretschmann scalar curvature has been shown for restricted initial data |136) . 
In the case of hyperbolic symmetry causal future geodesic completeness has been established by 
Rein |140j when the initial data are small. The plane and hyperbolic symmetric cases with a positive 
cosmological constant are analyzed in [185J. The authors show global existence to the future in 
areal time, and in particular they show that the spacetimes are future geodesically complete. The 
positivity of the cosmological constant is crucial for the latter result. A form of the cosmic no-hair 
conjecture is also obtained in [185 . It is shown that the de Sitter solution acts as a model for the 
dynamics of the solutions by proving that the generalized Kasner exponents tend to 1/3 as i — )• cxi, 
which in the plane case is the de Sitter solution. 

4.2.2 Gowdy and symmetric spacetimes 

The first study of spacetimes admitting a two-dimensional isometry group was carried out by 
Rendall |155j in the case of local symmetry. For a discussion of the possible topologies of these 
spacetimes we refer to the original paper. In the model case the spacetime is topologically of the 
form M X , and to simplify our discussion later on we write down the metric in areal coordinates 
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for this type of spacetime: 



g = e^^"-^) (-a + dO^) + e'^^t^ [dy + H d9 + M dtf 

+e^^[dx + Ady+{G + AH)de+{L + AM)dt]^. (52) 

Here the metric coefficients ij, U, a, A, H, L, and M depend on t and 6 and 6,x,y & S^. In [155] 
CMC coordinates are in fact considered rather than areal coordinates. Under the hypothesis that 
there exists at least one CMC hypersurface, Rendall proves for general initial data that the past 
of the given CMC hypersurface can be globally foliated by CMC hypersurfaces and that the mean 
curvature of these hypersurfaces blows up at the past singularity. The future direction was left 
open. The result in |155j holds for Vlasov matter and for matter described by a wave map. That 
the choice of matter model is important was shown in [154] , where a non-global existence result for 
dust is given, which leads to examples of spacetimes |101j that are not covered by a CMC foliation. 

There are several possible subcases to the symmetric class. The plane case, where the 
symmetry group is three-dimensional, is one subcase and the form of the metric in areal coordinates 



is obtained by letting A = G = H = L = AI = and U = log</2 in Equation (52 1. Another 



subcase, which still admits only two Killing fields (and which includes plane symmetry as a special 



case) , is Gowdy symmetry. It is obtained by letting G = H = L = M = Qm Equation ( 52 ) . In |6] 
Gowdy symmetric spacetimes with Vlasov matter are considered, and it is proven that the entire 
maximal globally hyperbolic spacetime can be foliated by constant areal time slices for general 
initial data. The areal coordinates are used in a direct way for showing global existence to the 
future, whereas the analysis for the past direction is carried out in conformal coordinates. These 
coordinates are not fixed to the geometry of spacetime and it is not clear that the entire past has 
been covered. A chain of geometrical arguments then shows that areal coordinates indeed cover the 
entire spacetime. The method in [6] was in turn inspired by the work |37 for vacuum spacetimes, 
where the idea of using conformal coordinates in the past direction was introduced. As pointed out 
in |2S], the result by Henkel [95, guarantees the existence of one CMC hypersurface in the Gowdy 
case and, together with the global areal foliation in [B], it follows that Gowdy spacetimes with 
Vlasov matter can be globally covered by CMC hypersurfaces as well. The more general case of 
symmetry was considered in [26] , where global CMC and areal time foliations were established 
for general initial data. In these results, the question whether or not the areal time coordinate 



takes values in (0, oo) or in {Rq, oo), Rq > 0, was left open. As we pointed out in Section 4.2.1 this 
issue was solved by Weaver [192] for symmetric spacetimes with the conclusion that Ro — 0, 
if the density function / is not identically zero initially. In the case of symmetric spacetimes, 
with a positive cosmological constant, Smulevici |174j has shown global existence in areal time 
with the property that t G (0, cxi). 

The issue of strong cosmic censorship for symmetric spacetimes has been studied by Dafer- 
mos and Rendall using the methods, which were developed in the surface symmetric case described 
above. In .66j strong cosmic censorship is shown under the same restriction on the initial data 
that was imposed in the surface symmetric case, which implies that the data have non-compact 
support in the momentum variable. Their result has been extended to the case with a positive 
cosmological constant by Smulevici |173j . 



4.3 Cosmological models with a scalar field 

The present cosmological observations indicate that the expansion of the universe is accelerating, 
and this has influenced theoretical studies in the field during the last decade. One way to produce 
models with accelerated expansion is to choose a positive cosmological constant. Another way is to 
include a non-linear scalar field among the matter fields, and in this section we review the results 
for the Einstein- Vlasov system, where a linear or non-linear scalar field have been included into 
the model. 
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Lee considers in |107) the case where a non-linear scalar field is coupled to Vlasov matter. The 
form of the energy momentum tensor then reads 



Tap = T^i-- + V„0V^0 - ( -V^0V^0 + ) 5a/3. (53) 



Here (j) is the scalar field and y is a potential, and the Bianchi identities lead to the following 
equation for the scalar field: 

V''V^(f>^ ¥'{(/)). (54) 

Under the assumption that is a non-negative function, global existence to the future is 
obtained, and if the potential is restricted to the form 

V{cl>) = Foe-^*, 

where < c < 4:^/^^, then future geodesic completeness is proven. 

In |187j the Einstein- Vlasov system with a linear scalar field is analyzed in the case of plane, 



spherical, and hyperbolic symmetry. Here, the potential V in Equations (53 1 and (54) is zero. A 
local existence theorem and a continuation criterion, involving bounds on derivatives of the scalar 
field in addition to a bound on the support of one of the moment variables, is proven. For the 
Einstein scalar field system, i.e., when f — 0, the continuation criterion is shown to be satisfied in 
the future direction, and global existence follows in that case. The work [186] extends the result 
in the plane and hyperbolic case to a global result in the future direction. In the plane case when 
/ = 0, the solutions are shown to be future geodesically complete. The past time direction is 
considered in [188j and global existence is proven. It is also shown that the singularity is crushing 
and that the Kretschmann scalar diverges uniformly as the singularity is approached. 



4.4 Stability of some cosmological models 

In standard cosmology, the universe is taken to be spatially homogeneous and isotropic. This 
is a strong assumption leading to severe restrictions of the possible geometries as well as of the 
topologies of the universe. Thus, it is natural to ask if small perturbations of an initial data set, 
which corresponds to an expanding model of the standard type, give rise to solutions that are 
similar globally to the future? 

In a recent work, Ringstrom |162| considers the Einstein- Vlasov system and he gives an affir- 
mative answer to the stability question for some of the standard cosmologies. 

The standard model of the universe is spatially homogeneous and isotropic, has flat spatial 
hypersurfaces of homogeneity, a positive cosmological constant and the matter content consists of 
a radiation fluid and dust. Hence, to investigate the question on stability it is natural to consider 
cosmological solutions with perfect fluid matter and a positive cosmological constant. However, 
as is shown by Ringstrom, the standard model can be well approximated by a solution of the 
Einstein- Vlasov system with a positive cosmological constant. Approximating dust with Vlasov 
matter is straightforward, whereas approximating a radiation fluid is not. By choosing the initial 
support of the distribution function suitably, Ringstrom shows that Vlasov matter can be made 
to mimic a radiation fluid for a prescribed amount of time; sooner or later the matter will behave 
like dust, but the time at which the approximation breaks down can be chosen to be large enough 
that the radiation is irrelevant to the future of that time in the standard picture. 

The main results in [162' are stability of expanding, spatially compact, spatially locally homo- 
geneous solutions to the Einstein- Vlasov system with a positive cosmological constant as well as a 
construction of solutions with arbitrary compact spatial topology. In other words, the assumption 
of almost spatial homogeneity and isotropy does not seem to impose a restriction on the allowed 
spatial topologies. 
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Let us mention here some related works although these do not concern the Einstein- Vlasov 
system. Ringstrom considers the case where the matter model is a non-linear scalar field in |163| 
and |164) . The background solutions, which Ringstrom perturb and which are shown to be stable, 
have accelerated expansion. In |163) the expansion is exponential and in [164] it is of power law 
type. The corresponding problem for a fluid has been treated in jl65j and |175j , and the Newtonian 
case is investigated in |109| and [10] for Vlasov and fluid matter respectively. 
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5 Stationary Asymptotically-Flat Solutions 



Equilibrium states in galactic dynamics can be described as static or stationary solutions of the 
Einstein- Vlasov system, or of the Vlasov-Poisson system in the Newtonian case. Here we consider 
the relativistic case and we refer to the excellent review paper |141| for the Newtonian case. First, 
we discuss spherically-symmetric solutions for which the structure is quite well understood. On 
the other hand, almost nothing is known about the stability of the spherically-symmetric static 
solutions of the Einstein-Vlasov system, which is in sharp contrast to the situation for the Vlasov- 
Poisson system. At the end of this section a recent result 1181 on axisymmetric static solutions will 
be presented. 

5.1 Existence of spherically-symmetric static solutions 

Let us assume that spacetime is static and spherically symmetric. Let the metric have the form 

where r > 0, 6* £ [0, tt], e [0, 27r]. As before, asymptotic flatness is expressed by the boundary 
conditions 

lim A(r) = lim /i(r) = 0, 

and a regular center requires 

A(0) = 0. 



Following the notation in Section 3.1 the time-independent Einstein-Vlasov system reads 



e-2^(2rA,. - 1) + 1 = STrr^p, (56) 
e-2^(2r^^ + 1) - 1 = 87rrV (57) 



Recall that there is an additional Equation (39) of second order, which contains the tangential 



pressure pr, but we leave it out since it follows from the equations above. The matter quantities 
are defined as before: 



p{x) = / ^l + \v\^J{x,v)dv, (58) 
- ^ (59) 



p{x) = / Z'^^LJ^^ f{x,v) — 

Jr3^ r J ^1 + \v\' 



The quantities 

E := p/'Mv/TTH^, L = l^n^p - (x ■ w)2 = |x X wp, 

are conserved along characteristics. E is the particle energy and L is the angular momentum 
squared. If we let 

fix,v)^^E,L), (60) 

for some function $, the Vlasov equation is automatically satisfied. 

A common assumption in the literature is to restrict the form of $ to 

'fiE,L)^<j>iE)iL-Loy+, (61) 

where / > —1/2, Lq > and — max{x, 0}. If we furthermore assume that 

^(E) ^ (E ^ Eo)l, k>-l, 
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for some positive constant Eq , then we obtain the polytropic ansatz. The case of isotropic pressure 
is obtained by letting I = so that $ only depends on E. 

In passing, we mention that for the Vlasov-Poisson system it has been shown |35j that every 
static spherically-symmetric solution must have the form / = <&(£', L). This is referred to as Jeans' 
theorem. It was an open question for some time whether or not this was also true for the Einstein- 
Vlasov system. This was settled in 1999 by Schaeffer [169] . who found solutions that do not have 
this particular form globally on phase space, and consequently, Jeans' theorem is not valid in the 
relativistic case. However, almost all results on static solutions are based on this ansatz. 

By inserting the ansatz f(x, v) — ^{E, L) in the matter quantities p and p, a non-linear system 
for A and /i is obtained, in which 

e-'^^{2r\r - 1) + 1 = 87rr2G$(r,/i), 
e^^^(2r^^ + 1) - 1 = 8TTr^H^{r,fi), 

where 



Ji Jo - 1 - L/r2 

r) ^OO ^r^(£^ — 1) 

= — / / $(e^«e,L)Ve2 - 1 - L/r^ dLde. 

^ Ji Jo 

Existence of solutions to this system was first proven in the case of isotropic pressure in jl44j, 
and extended to anisotropic pressure in [134] . The main difficulty is to show that the solutions 
have finite ADM mass and compact support. The argument in these works to obtain a solution 
of compact support is to perturb a steady state of the Vlasov-Poisson system, which is known 
to have compact support. Two different types of solutions are constructed, those with a regular 
centre [1441 1134j , and those with a Schwarzschild singularity in the centre [134j . 

In [145j . Rein and Rendall go beyond the polytropic ansatz and obtain steady states with 
compact support and finite mass under the assumption that $ satisfies 

$(£;, L) = c{E - Eq)\L^ + 0{{Ea - eY+'')L^ as E ^ Ea, 

where k > —1, I > —1/2, k + 1 + 1/2 > 0, k < I + 3/2. This result is obtained in a more direct way 
and is not based on the perturbation argument used in |1441 1134| . Their method is inspired by a 



work on stellar models by Makino |114j , in which he considers steady states of the Euler-Einstein 
system. In |145l there is also a discussion about steady states that appear in the astrophysics 
literature, and it is shown that their result applies to most of these steady states. An alternative 
method to obtain steady states with finite radius and finite mass, which is based on a dynamical 
system analysis, is given in \76\ . 



5.2 The structure of spherically-symmetric steady states 

All solutions described so far have the property that the support of p contains a ball about the 
centre. In [138] Rein showed that steady states also exist whose support is a finite, spherically- 
symmetric shell with a vacuum region in the center. In |S] it was shown that there are shell solu- 
tions, which have an arbitrarily thin thickness. A systematic study of the structure of spherically- 
symmetric static solutions was carried out mainly by numerical means in |22j and we now present 
the conclusions of this investigation. 

By prescribing the value /i(0), the equations can be solved, but the resulting solution will in 
general not satisfy the boundary condition /i(oo) = 0, but it will have some finite limit. It is then 
possible to shift both the cut-off energy Eq and the solution by this limit to obtain a solution. 
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which satisfies /Lt(oo) = 0. A convenient way to handle the problem that Eq and /i(0) cannot both 
be treated as free parameters is to use the ansatz 

fix, v)^{l- E/Eo)l{L - Lo)V, k>0, l> -1/2, k < 31 + 7/2, 

as in |22) . This gives an equation for /i, which can be rewritten in terms of the function 

Eo 

In this way the cut-off energy disappears as a free parameter of the problem and we thus have the 
four free parameters k,l,Lo and y(0). The structure of the static solutions obtained in [52] is as 
follows: 

If Lg = the energy density can be strictly positive or vanish at r = (depending on I) but 
it is always strictly positive sufficiently close to r = 0. Hence, the support of the matter is an 
interval [0, Ri] with i?i > 0, and we call such states ball configurations. If Lq > the support is 
in an interval [Rq, Ri], Rq > 0, and we call such steady states for shells. 

The value y{0) determines how compact or relativistic the steady state is, and the smaller values 
the more relativistic. For large values, recall y(0) < I/Eq, a pure shell or a pure ball configuration 
is obtained, cf. Figure [T] for a pure shell. Note that we depict the behavior of p but we remark that 
the pressure terms behave similarly but the amplitudes of p and pr can be very different, i.e., the 
steady states can be highly anisotropic. 

k=0, 1=1 .5, Ljj=0.2, y(0)=0.4 

2 
1.8 
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o. 1 
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Figure 1: A pure shell 

For moderate values of ?;(0) the solutions have a distinct inner peak and a tail-like outer peak, 
and by making ?;(0) smaller more peaks appear, cf. Figure [2] for the case of ball configurations. 

In the case of shells there is a similar structure but in this case the peaks can either be separated 
by vacuum regions or by thin atmospheric regions as in the case of ball configurations. An example 
with multi-peaks, where some of the peaks are separated by vacuum regions, is given in Figure [3j 

A different feature of the structure of static solutions is the issue of spirals. For a fixed ansatz of 
the density function /, there is a one-parameter family of static solutions, which are parameterized 



32 



k=2, 1=12, y(0)=0.05 
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Figure 3: Multi-peaks of a shell 
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by 2/(0)- A natural question to ask is how the ADM mass M and the radius of the support R 
change along such a family. By plotting for each j/(0) the resulting values for R and M a curve 
is obtained, which reflects how radius and mass are related along such a one-parameter family of 
steady states. This curve has a spiral form, cf. Figure |4j It is shown in that in the isotropic 
case, where I = Lq = Q the radius-mass curves always have a spiral form. 




0.15 



Figure 4: [R, M) spiral for fc = 0, 10.5, Lq = 0, 0.01 < y(0) < 0.99 



5.3 Buchdahl-type inequalities 

Another aspect of the structure of steady states investigated numerically in 22\ concerns the Buch- 
dahl inequality. If a steady state has support in [Rq, Ri], then the ADM mass M of the configura- 



tion is AI = m,{Ri), where the quasi local mass m,[r) is given by Equation (50 1 in Schwarzschild 
coordinates. 



In view of the Schwarzschild metric (51), Schwarzschild asked already in 1916 the question: 



How large can 2M/R possibly be? He gave the answer |170j 2M/R < 8/9 in the special case 
of the Schwarzschild interior solution, which has constant energy density and isotropic pressure. 
In 1959 Buchdahl |41] extended his result to isotropic solutions for which the energy density is 
non-increasing outwards and he showed that also in this case 



2M 



< . 

9 



(62) 



This is sometimes called the Buchdahl inequality. Let us remark that the Buchdahl inequality 
can obviously be written as M/R < 4/9, but since it is the quantity 2AI/R, which appears in the 
Schwarzschild metric (51), it is common to keep the form of Equation (62). A bound on 2M/R 
has an immediate observational consequence since it limits the possible gravitational red shift of a 
spherically-symmetric static object. 

The assumptions made by Buchdahl are very restrictive. In particular, the overwhelming 
number of the steady states of the Einstein- Vlasov system have neither an isotropic pressure nor 
a non- increasing energy density, but nevertheless 2M/R is always found to be less than 8/9 in 
the numerical study [22^. Also for other matter models the assumptions are not satisfying. As 
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pointed out by Guven and O Murchadha [HI] , neither of the Buchdahl assumptions hold in a simple 
soap bubble and they do not approximate any known topologically stable field configuration. In 
addition, there are also several astrophysical models of stars, which are anisotropic. Lemaitre |110| 
proposed a model of an anisotropic star already in 1933, and Binney and Tremaine ^38j explicitly 
allow for an anisotropy coefficient. Hence, it is an important question to investigate bounds on 
2M/R under less restrictive assumptions. 

In [ID] it is shown that for any static solution of the spherically-symmetric Einstein equation, 
not necessarily of the Einstein- Vlasov system, for which p > 0, and 



the following inequality holds 



p + 2pT< np, (63) 

Mr) (1 + 217)^-1 

r - (l + 2f7)2 ■ ' 

Moreover, the inequality is sharp and sharpness is obtained uniquely by an infinitely thin shell 
solution. Note in particular that for Vlasov matter Q. = 1 and that the right-hand side then 
equals 8/9 as in the Buchdahl inequality. An alternative proof was given in |103j and their method 



applies to a larger class of conditions on p,p and pT than the one given in Equation (63). On the 
other hand, the result in [103] is weaker than the result in ^lOJ in the sense that the latter method 
implies that the steady state that saturates the inequality is unique; it is an infinitely thin shell. 
The studies [TU[|103j are of general character and in particular it is not shown that solutions exist to 
the coupled Einstein-matter system, which can saturate the inequality. For instance, it is natural to 
ask if there are solutions of the Einstein- Vlasov system, which have 2m /r arbitrarily close to 8/9. 
This question is given an affirmative answer in [8], where in particular it is shown that arbitrarily 
thin shells exist, which are regular solutions of the spherically-symmetric Einstein- Vlasov system. 
Using the strategy in [9^ it follows that 

2m(r) 8 
T^^9' 

in the limit when the shells become infinitely thin. 

The question of finding an upper bound on 2M/R can be extended to charged objects and 
to the case with a positive cosmological constant. The spacetime outside a spherically-symmetric 
charged object is given by the Reissner-Nordstrom metric 

ds^ ^-{l-^+'^)dt^ + {l-^ + dr' + r\de^ + sin2 

where Q is the total charge of the object. The quantity 1 — ^ + ^ is zero when r± = M ± 
M 2 — Q^, and r± is called the inner and outer horizon respectively of a Reissner-Nordstrom 
black hole. A Buchdahl type inequality gives a lower bound of the area radius of a static object 
and this radius is thus often called the critical stability radius. It is shown in [11 that a spherically- 
symmetric static solution of the Einstein-Maxwell system for which p > 0, p+2pT < p, and Q < M 
satisfy 

Note, in particular, that the inequality holds for solutions of the Einstein- Vlasov-Maxwell system. 



since the conditions above are always satisfied in this case. This inequality (65) implies that 
the stability radius is outside the outer horizon of a Reissner-Nordstrom black hole. In [7S] the 
relevance of an inequality of this kind on aspects in black-hole physics is discussed. In contrast 
to the case without charge, the saturating solution is not unique. An infinitely thin shell solution 
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does saturate the inequality ( 65 ) , but numerical evidence is given in [16j that there is also another 
type of solution, which saturates the inequality for which the inner and outer horizon coincide. 

The study in |13| is concerned with the non-charged situation when a positive cosmological 
constant A is included. The following inequality is derived 

for solutions for which p > 0, p + 2pT < p, and < AR^ < 1. In this situation, the question 
of sharpness is essentially open. An infinitely thin shell solution does not generally saturate the 
inequality but does so in the two degenerate situations AR^ = and AR^ = 1. In the latter 
case there is a constant density solution, and the exterior spacetime is the Nariai solution, which 
saturates the inequality and the saturating solution is thus non-unique. In this case, the cosmo- 
logical horizon and the black hole horizon coincide, which is in analogy with the charged situation 
described above where the inner and outer horizons coincide when uniqueness is likely lost. 



5.4 Stability 

An important problem is the question of the stability of spherically-symmetric steady states. At 
present, there are almost no theoretical results on the stability of the steady states of the Einstein- 
Vlasov system. Wolansky |195j has applied the energy-Casimir method and obtained some insights, 
but the theory is much less developed than in the Vlasov-Poisson case and the stability problem 
is essentially open. The situation is very different for the Vlasov-Poisson system, and we refer 
to |141) for a review on the results in this case. 

However, there are numerical studies [3T1 llOOl I171j on the stability of spherically-symmetric 
steady states for the Einstein- Vlasov system. The latter two studies concern isotropic steady states, 
whereas the first, in addition, treats anisotropic steady states. Here we present the conclusions 
of [H], emphasizing that these agree with the conclusions in |171l 1100] for isotropic states. 

To allow for trapped surfaces, maximal-areal coordinates are used, i.e., the metric is written in 
the following form in [2T] 

ds^ ^ -{a^ - a^p'^)dt^ + 2a'^pdtdr + a'^dr'^ + {dO'^ + sin^ 9 dcj)^) . 

Here the metric coefficients a, (3, and a depend on t € K and r > 0, a and a are positive, and the 
polar angles S [0, tt] and G [0, 27r] parameterize the unit sphere. Thus, the radial coordinate r is 
the area radius. A maximal gauge condition is then imposed, which means that each hypersurface 
of constant t has vanishing mean curvature. The boundary conditions, which guarantee asymptotic 
flatness and a regular center, are given by 

a(t,0) = a(t,oo) = a{t, oo) = 1. (66) 

Steady states are numerically constructed, and these are then perturbed in order to investigate 
the stability. More precisely, to construct the steady states the polytropic ansatz is used, cf. 
Section |5.1[ 

/(r, w, L) = $(£;, L) = (Eo - E)X{L - io)V- (67) 

By specifying values on Eq, Lq and a{0) steady states are numerically constructed. The distribu- 
tion function fs of the steady state is then multiplied by an amplitude A, so that a new, perturbed 
distribution function is obtained. This is then used as initial datum in the evolution code. We 
remark that also other types of perturbations are analyzed in [21. . 

For k and I fixed each steady state is characterized by its central red shift Zc and its fractional 
binding energy Ei,, which are defined by 

1 Cf, 

= —77^ - 1, Eb = — , where = Afo - M. 
a[0) Mo 
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Here 

pOO pOO pOO 

Mo=47rM / / a{t,r)f{t,r,w,L)dLdwdr 

Jo J-oo Jo 

is the total number of particles, which, since all particles have rest mass one, equals the rest mass 
of the system. M is the ADM mass given by 



M 



^47rp(t, r) + -^li'^it, r) ) r^dr, 



where At = fi/ra. Both Mq and M are conserved quantities. The central redshift is the redshift of a 
photon emitted from the center and received at infinity, and the binding energy eb is the difference 
of the rest mass and the ADM mass. In Figure [5] and Figure |6] the relation between the fractional 
binding energy and the central redshift is given for two different cases. 
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Figure 5: fc = 0, / = 0, 0.1 

The relevance of these concepts for the stability properties of steady states was first discussed by 
Zel'dovich and Podurets |197| . who argued that it should be possible to diagnose the stability from 
binding energy considerations. Zel'dovich and Novikov |196| then conjectured that the binding 
energy maximum along a steady state sequence signals the onset of instability. 

The picture that arises from the simulations in [2T] is summarized in Table [l] Varying the 
parameters k,l and Lq give rise to essentially the same tables, cf. 21j. 

If we first consider perturbations with A > 1, it is found that steady states with small values 
on Zc (less than approximately 0.40 in this case) are stable, i.e., the perturbed solutions stay in a 
neighbourhood of the static solution. A careful investigation of the perturbed solutions indicates 
that they oscillate in a periodic way. For larger values of Zc the evolution leads to the formation 
of trapped surfaces and collapse to black holes. Hence, for perturbations with A> 1 the value of 
Zc alone seems to determine the stability features of the steady states. Plotting Ei, versus Zc with 
higher resolution, cf. |21| . gives support to the conjecture by Novikov and Zel'dovich mentioned 
above that the maximum of Ei, along a sequence of steady states signals the onset of instability. 

The situation is quite different for perturbations with A < 1. The crucial quantity in this 
case is the fractional binding energy E^. Consider a steady state with Eb> Q and a perturbation 
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Figure 6: Zc 0, ^ 3/2, Lq = 0.1 

with A < 1 but close to 1 so that the fractional binding energy remains positive. The perturbed 
solution then drifts outwards, turns back and reimplodes, and comes close to its initial state, and 
then continues to expand and reimplode and thus oscillates, cf. Figure [7j 

In |171j it is stated (without proof) that if > the solution must ultimately reimplode and 
the simulations in [5T] support that it is true. For negative values of Ei,, the solutions with A < 1 
disperse to infinity. 

A simple analytic argument is given in |21j . which relates the question, whether a solution 
disperses or not. It is shown that if a shell solution has an expanding vacuum region of radius 
R{t) at the center with R{t) — > oo for f — )■ oo, i.e., the solution disperses in a strong sense, then 
necessarily Mq < M, i.e., Ef, < 0. 

5.5 Existence of axisymmetric static solutions 

As we have seen above, a broad variety of static solutions of the Einstein- Vlasov system has 
been established, all of which share the restriction that they are spherically symmetric. The 
recent investigation [TS] removes this restriction and proves the existence of static solutions of the 
Einstein- Vlasov system, which are axially symmetric but not spherically symmetric. From the 
applications point of view this symmetry assumption is more "realistic" than spherical symmetry, 
and from the mathematics point of view the complexity of the Einstein field equations increases 
drastically if one gives up spherical symmetry. Before discussing this result, let us mention that 
similar results have been obtained for two other matter models. In the case of a perfect fluid, Heilig 
showed the existence of axisymmetric stationary solutions in [92 . These solutions have non-zero 
angular momentum since static solutions arc necessarily spherically symmetric. In this respect the 
situation for elastic matter is more similar to Vlasov matter. The existence of static axisymmetric 
solutions of elastic matter, which are not spherically symmetric, was proven in [1]. Stationary 
solutions with rotation were then established in [2j. 

Let us now briefly discuss the method of proof in [18] , which relies on an application of the 
implicit function theorem. Also, the proofs in [921 [Tl make use of the implicit function theorem. 
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Figure 7: = 0.47, = 0.04, A = 0.99, T = 90.0 
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Table 1: fc = and Z = 1/2. 
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unstable 


unstable 



but apart from this fact the methods are quite different. 

The set-up of the problem in [18' follows the work of Bardeen [STl , where the metric is written 
in the form 

for functions 1^,3,(1 depending on p and z, where t G M, p e [0, cxd[, z G M and ip € [0, 27r]. The 
Killing vector fields dt and d(j, correspond to the stationarity and axial symmetry of the spacetime. 
Solutions are obtained by perturbing off spherically symmetric steady states of the Vlasov-Poisson 
system via the implicit function theorem and the reason for writing iz/c^ in the metric, instead 
of V, is that v converges to the Newtonian potential Un of the steady state in the limit c — >■ (30. 
Asymptotic flatness is expressed by the boundary conditions 



lim v{p^z) = lim fi(p,z) = 0, lim B{p,z) = 1. 

|(p,z)|-i.oo \(p,z)\^oo |(p,z)|-i.oo 



(69) 



In addition the solutions are required to be locally flat at the axis of symmetry, which implies the 
condition 

^(0,z)/c2 +^(0,z) = lnB(0,z), zeR. (70) 

Let us now recall from Section |5.1| the strategy to construct static solutions by using an ansatz of 
the form 

where E and L are conserved quantities along characteristics. Due to the symmetries of the 



metric (68 1 the following quantities are constant along geodesies: 
E 



'9{d/dt,p-) 



+ C-2 (e2f (pl)2 + e2M(p2)2 _,,p2^2g-2i./c2(p3)2)^ 



L := gid/d^.p") = p" B^e-^''/^'/- 



(71) 
(72) 



Here p° are the canonical momenta. E can be thought of as a local or particle energy and L is the 
angular momentum of a particle with respect to the axis of symmetry. For a sufficiently regular 
4> the ansatz function / satisfies the Vlasov equation and upon insertion of this ansatz into the 



definition of the energy momentum tensor (32 ) the latter becomes a functional T^p — Tap{v, B, p) 
of the unknown metric functions v, B, p. It then remains to solve the Einstein equations with this 
energy momentum tensor as right-hand side. The Newtonian limit of the Einstein- Vlasov system is 
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the Vlasov-Poisson system and the strategy in |18j is to perturb off sphericahy symmetric steady 
states of the Vlasov-Poisson system via the impHcit function theorem to obtain axisymmetric 
solutions. Indeed, the main result of |18| specifies conditions on the ansatz function $ such that 
a two parameter (7 and A) family of axially-symmetric solutions of the Einstein- Vlasov system 
passes through the corresponding spherically symmetric, Newtonian steady state, whose ansatz 
function we denote by 4>. The parameter 7 = 1 /c^ turns on general relativity and the parameter A 
turns on the dependence on L. Since L is not invariant under arbitrary rotations about the origin 
the solution is not spherically symmetric if / depends on i. It should also be mentioned that 
although 7 is a priori small, which means that c is large, the scaling symmetry of the Einstein- 
Vlasov system can be used to obtain solutions corresponding to the physically correct value of c. 
The most striking condition on the ansatz function <&, or rather on the ansatz function (p of the 
corresponding Vlasov-Poisson system, needed to carry out the proof is that it must satisfy 



An important argument in the proof is indeed to justify that there are steady states of the Vlasov- 
Poisson system satisfying this condition. 

It is of course desirable to extend the result in [TH] to stationary solutions with rotation. More- 
over, the deviation from spherically symmetry of the solutions in jl8j is small and an interesting 
open question is the existence of disk-like models of galaxies. In the Vlasov-Poisson case this has 
been shown in |74) . 
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